Fluctuations of the energy levels number at the mobility edge. 
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We present a detailed analysis of the behavior of two-level 
correlation function R(s) in the disordered sample. We show 
that in the vicinity of the Anderson transition as well as in 
2d, the variance of the number of levels in an energy band of 
a width E has a linear behavior as a function of E. This is 
related with an "anomalous" contribution to the sum rule for 
R(s), which comes from the "ballistic" region of s. 
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The problem of energy level correlation in quantum 
systems has been studied intensively after the pioneer 
work of Wigner |Q. The random matrix theory devel- 
oped by Wigner and Dyson B is known to describe the 
statistical properties of spectra of various complex sys- 
tems, such as nuclei, atoms or molecules. Gor'kov and 
Eliashberg || put forward an assumption that the ran- 
dom matrix theory can be also applied to the problem of 
energy level statistics of a quantum particle subject to a 
random potential. This hypothesis was proven by Efetov 
101 , who studied the two-level correlation function 



(i) 



where v{e) is the density of states at the energy e, (. . .) de- 
note averaging over realizations of the disorder, s = to/ A., 
and A = {{v)V)~ 1 is the mean level spacing. As was 
shown in Ref. 0, the correlator R(s) for a disordered 
mesoscopic metallic sample is described by the Wigner- 
Dyson distribution for u> <C E c , where E c is the Thouless 
energy; E c — D/L 2 , with D being the diffusion con- 
stant and L the sample size. Corrections to the Wigner- 
Dyson distribution in this region were studied in a recent 
paper B. These corrections are of order 1/g 2 , where 
g = Ec/A 1 is the dimensionless conductance of the 
sample. For lu ^> E c , the behavior of the correlation func- 
tion changes because the corresponding time scale cj _1 is 
smaller than the time E~ l the particle needs to diffuse 
through the sample. The smooth (non-oscillating) part 
of the correlation function in this region was calculated 
in H by means of the diffusion perturbation theory. 

All these results were obtained for the case of weak 
disorder (clean enough sample), when j > 1. As is well 
known, when the degree of disorder increases, the Ander- 
son localization phenomenon takes place. In 3d, there is 
a critical value of conductance g = g* ~ 1, which cor- 
responds to the metal insulator transition. The energy 
level statistics near this mobility edge is different from 



that of a good metal and has attracted a large inter- 
est recently. The first study of this problem was under- 
taken by Altshuler et al. , who considered the variance 
{SN 2 (E)) = (N 2 (E)) - (N(E)) 2 of the number of levels 
within a band of a width E. This quantity is related to 
the two-level correlator (Q) via 

f (N(E)) 

(SN(E) 2 ) = / ((N(E))) - \s\)R(s)ds , (2) 

J-(N{E)) 

or, equivalently, 

— — — (SN(E) 2 ) = / R( S )ds , (3) 

d(N(E)) J-{N(E)) 

where (N(E)) = E/A. As was shown by Dyson 0), the 
random matrix theory leads to the logarithmic behavior 
of the variance (8N 2 ) WD ~ ln(A) for (A) > 1. In the 
opposite situation, when all energy levels are completely 
uncorrelated (known as the Poisson statistics), one gets 
{5N 2 )p — (N). Based on some scaling considerations as 
well as on numerical simulations, Altshuler et al. 
forward a conjecture that at the critical point 



Put 



(SN 2 ) ~ k{N) 



(4) 



where < k < 1 is certain numerical coefficient. More 
recently, Shklovskii et al. ||] introduced the conception 
of new universal statistics at the mobility edge. In Ref. 
|1 the correlator R(s) at the mobility edge was studied 
by means of perturbation theory combined with scaling 
assumptions about a form of the diffusion propagator. It 
was found that for s 3> 1, 



R(s) oc s- 2+7 



(5) 



where 7 < 1 is certain critical index. Using this result, 
the relation (0) and the sum rule 



R(s)ds = 



(G) 



(which follows from the conservation of the number of 
energy levels), the authors of U concluded that 



(SN 2 ) oc N"< 



(7) 



with 7 < 1, in contradiction with Ref. 0]. 

In this paper, we present a thorough analysis of the 
behavior of the correlator R(s) in all range of values of s 
and for different strength of the disorder, corresponding 
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to a good metal, a critical region and a critical point. We 
show that the asymptotic behavior (||) of the correlator 
R(s) at s ^> 1 does not imply the absence of the linear 
term (Q). It was assumed in || that the universal part 
of the correlator R(s) (which is the one surviving in the 
limit E/A = (N) = const, L — > oo) satisfies the sum rule 
(^J). We will show however that the sum rule is fulfilled 
only if all contributions are taken into account, including 
the non-universal contribution of the "ballistic" region 
lu ~ 1/t, where r is elastic mean free time. For this 
purpose, we will estimate the contributions to the sum 
rule from all regions of the variable s. 

In fact, for the conventional model of a particle in a 
random potential the definition (jl]) of the correlator and 
the sum rule relation (JE)) should be modified, when the 
vicinity of critical point is considered. The reason is that 
the Anderson transition point corresponds to a strong 
disorder regime ct ~ 1, so that the condition lu <~ 1/t 
implies lu <~ e. On the other hand, the density of states 
v{e) can be considered as a constant only for small varia- 
tions of energy lu <C e. This means that variation of v (E) 
should be taken into account in (|l]), (Q). Besides, the con- 
dition er ~ 1 leads to a breakdown of the perturbation 
theory, that complicates the analysis of the "ballistic" 
region contribution. To get rid of these problems, we 
consider here a different microscopic model which has ex- 
actly the same universal part of R(s), but whose density 
of states does not change within the range of lu ~ 1 jr. 
This is so-called n-orbital Wegner model |l(J , which can 
be considered as a system of metallic granulae forming a 
d-dimensional lattice, each granula being coupled to its 
nearest neighbours. In the limit n 3> 1 this model can 
be mapped onto a supersymmetric er-model defined on a 
lattice. The action of this er-model reads lfOl.1121 



S{Q} = ~ E S^QiQj + e E Str ^ A 

(ij) ' 



(8) 



where Qt = T+ATf are 4x4 supermatrices defined on 
sites i of a d-dimensional lattice with a lattice spac- 
ing a, Ti belongs to the coset space f/(l, 1|2), A = 
diag{l, 1,-1,-1} and Str denotes the supertrace. Sum- 
mation in the first term of eq.(]^) goes over pairs (ij) of 
nearest neighbours. The parameters 7 and e are related 
to the classical diffusion constant D, the density of states 
v and the frequency lu as follows: 



7 = TTfDa' 



d-2 



e = — iunua /2 



(9) 



To find a detailed description of the model and relevant 
mathematical objects, a reader is referred to Refs. P,|TT|. 

The two-level correlator can be expressed through a 
correlation function of this er-model in the following way: 

R{s) = Re J JJ DQ A StrQM 

j [Li j J 

x exp(-S{Q}) , (10) 



where k = diag{l, — 1, 1, — 1} and V is the system vol- 
ume. An important feature of this lattice er-model is 
that the corresponding correlator (10) satisfies exactly 
the sum rule ([|) (r|. As will become clear below, the 
continuum version of the er-model does not possess this 
property: there is a deficiency of the sum rule for it which 
is related with the contribution of the domain of lu close 
to an ultraviolet cut-off. This phenomenon is similar to 
what is known as quantum anomalies in the quantum 
field theory (see e.g. ]1^| ). To prove the sum rule for 
the correlator ( JlC| ) on a lattice, we substitute eq.(|8|) into 
cq.dTcl) and integrate over lu. This gives 



J R(s)ds = Re J Y[DQj j ^ ^ StrQ.fcA 
x ex P j ~\ E S^QiQj j E St ^A) 



- 1 



(11) 



In the conventional way |4|,|1 lj , the Q-matrices can be 
parametrized by 2 "eigenvalues" 1 < Ai < 00, — 1 < 
A2 < 1 and a certain set of "angular" variables. We do 
not need an explicit form of this parametrization here, 
but only the fact that StrQA = 2(Ai - A 2 ) > 0, so that 
StrQA = only when Ai = A2 = 1. Therefore, the 



hyperplane 
function in (I 




Ag ) = determined by the 5- 



only at one point: X\ l> = X\ l> = 1 for all i, 



LID touches the actual domain of integration 

\^ = Ao = 1 f° r au h that leads to 
vanishing of the integral |15[ . 

Let us stress that in the region lu <C D/a 2 = E c (a) 
the correlator R(s) is universal, i.e. does not depend on 
microscopical details of the model. The region lu ~ D /a 2 
plays a role analogous to the ballistic one lu ~ 1/t in the 
case of a usual model of particle in a random potential. 
Despite the non-universality of the correlation function 
R nu (s) in this "ballistic" domain, the corresponding inte- 
gral contribution /„„ to the sum rule is universal, because 
it determines, according to eq.(^), the sum rule deficiency 
for the universal part R u (s): 

~T" Ryu , 

lu = J R u (s)ds ; I nu = J R nu (s)ds (12) 

As will be seen below, when the system is close to the An- 
derson transition the two regions of the variable s domi- 
nating the integrals I u and I nu respectively are separated 
by a parametrically broad domain of s which gives a neg- 
ligible contribution to the sum rule. For energy band 
width E lying in this region, we get from eqs.(0), (|12|) 



(5N{Ef) ~ I U {N{E)} = -I nu (N{E)) 



(13) 



i.e. just the linear term (|J), which was conjectured in 
Ref. |Q. In Rcf. j| the "anomalous" contribution /„„ 
was lost, that led to a conclusion about the absence of 
this linear term. 
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We come now to the analysis of the correlator R(s) 
and of the sum rule in various situation. We start from 
the case of a good metal, then we consider the critical 
point and the critical region cases. 

1. Good metal. Here the following three regions with 
different behavior of R(s) can be found: 

A. "Wigner-Dyson" (WD) region: lu <C E c . The cor- 
relator R(s) in this region is given by the random matrix 
theory result Rw d (s) plus a, 1/g 2 correction B : 



R(s) ~ Rwd(s) 



a d -2/ \ 

„ „ sin (tts) 



(14) 



where ad > is certain numerical coefficient. Corre- 
sponding contribution to the sum rule can be estimated 
as 



Ia 



R(s) oc +l/g 



(15) 



In all estimates we omit numerical factors of order unity, 
keeping only a sign and a parametric dependence. 

B. "Altshuler-Shklovskii" (AS) region: E c < lu < 
E c (a). Here E c {a) = D/a 2 is the Thouless energy at 
the scale of lattice spacing a, which plays a role of the 
ultraviolet cut-off for the diffusion theory. We have M : 



R(s) oc +g 



-d/2 s d/2-2 



(16) 



and consequently, 



D/a 2 A 



R(s)ds cx g- d/2 



D 
~a^A 



oc 



-l/5(o) 



(17) 



where g(a) = g(a/L) d ~ 2 is the conductance at the scale 
a. Note that for the case of a particle in a random poten- 
tial, the following substitutions should be done: a — » I; 
E c (a) -> E c (l) = 1/r; l/g{a) -> l/g(l) cx (er)^ 1 . 

C. "Ballistic" region: ui > E c (a). (As was already 
noted, this region is not really ballistic in the lattice a- 
model, but plays exactly the same role as the ballistic 
one in the usual case.) To find the correlator R(s) in this 
region, we can neglect to the first approximation the first 
term in eq.(0), that gives 



L\ d 1 



I c - / i?(s)ds oc -l/5(a) 



(18) 
(19) 



The contribution (g_5J) to the sum rule is dominated by 
the region lu ~ A, whereas the contributions (p7|), jTg| ) 
are dominated by a vicinity of the ultraviolet cut-off 
lu ~ D/a 2 . In d = 3 3 > 5(1)7 and these non-universal 
contributions > are much larger (by absolute value) 
than I a, but they should cancel each other according 
to the sum rule. The situation is more interesting in 



d = 2, where the contribution ([l6]), flT7| ) of the AS region 
is absent in view of special analytical properties of the 
diffusion propagator In this case g(a) = g, so that 
the contributions Ia and Ic are of the same order, in 
agreement with the sum rule which prescribes their sum 
to be zero. We can identify then Ia with I u and Ic with 
I nu in eq.(|l2|), so that I u = —I nu oc I/5. According to 
eq.(|13), this leads to a linear term for the variance: 



(5N 2 ) cx (l/g)(N) 



E C ^E^. D/a 2 



(20) 



but the corresponding coefficient is of order 1 / g <C 1 . 

With an increase in disorder strength, the coupling 
constant 7 in eq.(||) decreases. When it approaches a 
critical value 7 C , corresponding to the metal-insulator 
transition, the correlation length £ becomes large: £ 3> a. 
Depending on the relation between £ and the system size 
L, the following two situations may take place. 

2. Metal in critical region: a <C £ <C L. In this case 
g(a) ~ but g = g{L) ^> 5*, where 3* is the value of 
conductance in the critical point |^9| . We find then as 
much as four different regions of s for the correlator R(s): 

A. WD region: to <C E c . The correlator R{s) and the 
contribution to the sum rule for this region are given by 
the same eqs.(14), (|l5|), as for a good metal. 

B. AS region: E c < u < E c (£) = 5*A 5 , where E c (£) 
and are Thouless energy and level spacing for a piece 
of the sample with linear size £. The correlator i?(s) is 
given by eq. ([l6|) , and we find 



S»A £ /A 



R(s)ds cx g-^ 2 ,?, I -r 5 - ) oc +l/.g 



A c \ 



(21) 

where we have used the relations A^/A = (L/£) d and 
.9/3, cx (L/£) d - 2 . 

C. "Kravtsov-Lerner-Altshuler-Aronov" (KLAA) re- 
gion: E c (£) = 3* A j <C u <C ff,A,, e E c (a). The correla- 
tor i?(s) in this region was studied in Ref. (l7|; the result 
reads: 



AGO c - 

and consequently 



-.9* 



(A^/A) 1 - 



-2+7 



rg.A a /A 

Ic — R(s)ds oc —1/g* 

Jg*A i /A 



(22) 



(23) 



D. "Ballistic" region: cj > E c (a). Here the equations 
©, © hold, yielding 



7_d cx -I/5* 



(24) 



The contribution I a is dominated by the domain to ~ 
A, the contributions Ib, Ic by w ~ A^ and finally the 
contribution Iu comes from the "non-universal" region 
lu ~ E c (a). Therefore, this last contribution determines 



i nu in eq.(12), and according to eq.(13) we find the linear 
behavior of the variance 
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(6N 2 ) cx (l/g»){N) ; E (£) < E « £ c (a) (25) 

with the coefficient which is of order unity in 3d. 

3. Finally, a system is at the critical point, when 
£ 3> £. In this case, = g(a) = g*. The AS region 
disappears, and we have, in full analogy with previous 
estimates: 

A. WD region: uj < E c = A. 

I A oc +1/.9, (26) 

B. KLAA region: g*A < cj <C A a g*. 

R(s) cx -^s" 2 ^ ; Ib oc (27) 

C. "Ballistic" region, co > A a g*. 

/c oc -l/<?* (28) 

Again the same conclusion, as for the critical region 
case, can be drawn: the "ballistic" contribution ( |28| ) can 
be identified as in eq. (U), that gives: 

(5N 2 ) cx (l/g*)(N) ; £ c < £ < £ c (a) (29) 

This concludes our analysis. In fact, all obtained re- 
sults can be summarized in the following relation: 

/u / A p 
R{s)ds = - R(s)ds oc l/g(L u ) , (30) 

-u/A J\s\>uj/A 

where is defined by the condition E C (L^) = ui, and 
g{LiS) is the conductance of a piece of the sample of the 
size L u . Equation ( |30| ) holds for all to such that L u < L. 
Its physical meaning is quite transparent: a portion of a 
state "expelled" from the volume with linear size L u is 
inverse proportional to its conductance. 

In conclusion, we have demonstrated that the levels 
number variance (||) shows a linear behavior ((2^), J29|) 
in the vicinity of Anderson transition, in agreement with 
an earlier prediction fif]. The corresponding coefficient 
is proportional to 1/g*, with g» being the critical value 
of conductance. We have shown that this linear term is 
related to an "anomalous" contribution to the sum rule 
for two-level correlation function, which originates from 
"ballistic" region of frequencies providing an ultraviolet 
cut-off for the diffusion theory 

To avoid possible confusion, let us note that one should 
distinguish the asymptotic behavior of (SN 2 ) from that 
of another frequently studied quantity - spacing distri- 
bution P(s). The former is by definition the variance of 
the distribution function of number of levels in an energy 
band of given width, whereas the latter is determined by 
the far tail of this distribution function. The asymptotic 
behavior of P(s) at the mobility edge was found in Ref. 
Hi to be P(s) ~ exp(-s 2 -T). We believe that the va- 
lidity of this result is not affected by the present study. 

We are grateful to V.E.Kravtsov for numerous discus- 
sions which stimulated appearence of this work. A.G.A. 



thanks B.Shapiro for attracting his attention to the prob- 
lem. A.D.M. acknowledges discussions of the results with 
Y.Gefen, F.Izrailev and Y.V.Fyodorov and is grateful to 
the International Centre for Theoretical Physics in Tri- 
este and the Institute for Nuclear Theory at the Univer- 
sity of Washington for the kind hospitality and to the 
Alexander von Humboldt Foundation for support. 



' Permanent address: Petersburg Nuclear Physics Insti- 
tute, 188350 Gatchina, St. Petersburg, Russia. 
[1] E.P.Wigner, Proc. Cambridge Philos. Soc. 47, 790 
(1951). 

[2] F.J.Dyson, J.Math.Phys. 3, 140 (1962). 

[3] L.P.Gor'kov and G.M.Eliashberg, Zh.Eksp.Theo-r.Fiz., 

48, 1407 (1965) [Sov. Phys. JETP, 21, 940 (1965)]. 
[4] K.B.Efetov, Adv. Phys., 32, 53 (1983). 
[5] V.E.Kravtsov, A.D.Mirlin, submitted to Phys.Rev.Lett. 
[6] B.L.Altshuler, B.I.Shklovskii, Zh.Eksp.Theor.Fiz. 91, 

220 (1986) [Sov.Phys. JETP 64, 127 (1986)]. 
[7] B.L.Altshuler, I.K.Zharekeshev, S.A.Kotochigova and 
B.I.Shklovskii, Zh.Eksp.Theor.Fiz. 94, 343 (1988) 
[Sov.Phys. JETP 67, 625 (1988)]. 
[8] B.I.Shklovskii, B.Shapiro, B.R. Sears, P.Lambrianides 

and H.B. Shore, Phys.Rev. B 47, 11487 (1993). 
[9] V.E.Kravtsov, I.V.Lerner, B.L.Altshuler, A.G.Aronov, 
Phys. Rev. Lett. 72, 888 (1994). 
[10] F.Wegner, Phys.Rep. 67, 15 (1980). 
[11] M.R.Zirnbauer, Nucl. Phys. B 265, 375 (1986). 
[12] In fact, Eqs.(|)-(0) are written for the case of cr-model 
of unitary symmetry, corresponding to broken time rever- 
sal invariance; all results are however completely analo- 
gous for orthogonal and symplectic symmetry cases. 
[13] The idea that the sum rule can be proven directly for the 

o"-model was suggested by V.E.Kravtsov. 
[14] C.Itzykson and J.-B.Zuber, Quantum Field Theory, 

McGraw-Hill, 1980. 
[15] In fact, the situation is slightly more complicated because 
of boundary terms which are known to appear when one 
calculates an integral over the super coset space [jll| . An 
integral over each matrix Qi will be given by a sum of 
such boundary term and a regular integral over the eigen- 
values Aj l ',Aj . The contribution to eq.(p4|) originating 
from boundary terms for all Qi will then cancel the dis- 
connected part of the correlator ( — 1 in brackets in (|ll|)). 
To all other terms the argumentation presented above can 
be applied. This proves vanishing of the integral ([n]). 
[16] A. Altland, Y.Gefen, Phys.Rev.Lett. 71, 3339 (1993); 

V.E.Kravtsov and I.V.Lerner, unpublished. 
[17] A.G.Aronov, V.E.Kravtsov, I.V.Lerner, submitted to 
Phys.Rev.Lett. 

[18] A.G.Aronov, V.E.Kravtsov, I.V.Lerner, JETP Lett. 59, 
40 (1994). 

[19] In 3d the value of g* is of course of order unity. However, 
if one consider formally d = 2 + e with e <C 1, then 
g» ~ 1/e and is therefore parametrically large. That is 
why it is reasonable to keep g, as a parameter in all 
estimates. 



4 



